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Abstract 

The dissociation of molecules, even the most simple hydrogen molecule, cannot be described 
accurately within density functional theory because none of the currently available functionals 
accounts for strong on-site correlation. This problem has led to a discussion of properties that 
the local Kohn-Sham potential has to satisfy in order to correctly describe strongly correlated 
systems. We derive an analytic expression for this potential at the dissociation limit and show that 
the numerical calculations for a one-dimensional two electron model system indeed approach and 
reach this limit. It is shown that the functional form of the potential is universal, i.e. independent 
of the details of the system. 

PACS numbers: 31.15.E-,31.15.V-,31.15.vn 
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I. INTRODUCTION 



Over the years the improvement in exchange- correlation (xc) functional has made density 
functional theory (DFT)[1, 2] the tool of choice to accurately study and predict properties 
of many-electron systems. Applications range from atoms to molecules and nanostructures, 
biomolecules and solids and cover diverse topics such as theoretical spectroscopy, e.g. op- 
tical, energy loss, and time-resolved spectroscopy, electron transport, light induced phase 
transitions, photochemistry, and electrochemistry [3, 4]. Despite this success major basic 
challenges remain that usually are manifestations of strong, static and dynamic, electron 
correlations [5]. Van der Waals interactions, the localization in strongly-correlated systems, 
open-shell molecules, and molecular dissociation are poorly accounted for by present func- 
tional [5, 6]. A general measure of inter-electron correlations is the ratio of the kinetic 
energy to the potential energy of the Coulomb interaction between electrons. While the 
kinetic energy is lowered by delocalization of electrons over the system, the Coulomb repul- 
sion works in the opposite direction trying to keep electrons far from each other and thus 
favoring the tendency to localization. In the condensed matter context this interplay of two 
opposite tendencies is commonly pictured in terms of the Hubbard on-site correlations that 
suppress tunneling of particles between atoms and lead to localization of electrons on lattice 
sites (or groups of sites). Strong Hubbard correlations are responsible for the dissociation of 
molecules, the physics of Mott insulators, non-itinerant magnetism in most of the magnetic 
dielectrics, the Coulomb blockade in quantum transport, etc. The failure of the common 
DFT-functionals to capture the effects of Hubbard correlations led to the development of the 
LDA-I-U method [7] and its more elaborated counterpart, the dynamical mean- field theory 
(DMFT) [8], to describe strongly correlated systems. On the other hand, it is absolutely 
clear that DFT being an "in principle exact" theory should be capable to describe the regime 
of strong correlations provided the proper xc potential is known. In this realm, it is funda- 
mental to increase the knowledge of relations, fulfilled by the exact xc potential, in order to 
move forward on the road towards the ultimate functional, the "holy-grail of DFT" . 

In the present work, we consider a prototypical example of a physical behavior governed 
by strong Hubbard correlations - the dissociation of diatomic molecules, and discuss exact 
features of the xc potential v^c necessary to describe the correlation-driven electron localiza- 
tion happening in the dissociation limit. One such feature is well known - in the dissociation 
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of heteroatomic molecules the Kohn-Sham (KS) potential Vg acquires a step in between the 
fragments to adjust the ionization potentials [9]. The value of this step is universal and is 
simply given by the difference of the ionization potentials of the two fragments of the disso- 
ciated molecule. Apparently, the presence of this step is necessary to prevent an unphysical 
fall of electrons to the fragment with a higher ionization potential. However, as we show 
below, it is not sufficient to correctly describe the dissociation, i.e. the strongly correlated, 
limit. In fact, in this limit the xc potential acquires a nontrivial structure even for the most 
simple homoatomic molecules, such as H2. 

An important step in understanding the behavior of the xc potential in the dissociation 
limit has been made in a series of works by Baerends and co-authors [10-14] who recon- 
structed the xc potential of a number of stretched diatomic molecules from an accurate 
many-body configuration interaction (CI) ground-state wave function. They noticed numer- 
ically that, in addition to the step, v^c also shows a peak structure around the middle point 
between the two atoms [10-13]. A subsequent analysis has shown that the peak in v^c is 
probably a general feature of the dissociation limit, which contradicts the common LCAO 
form of the molecular orbital, but can be reasonably well reproduced assuming that the 
two-electron wave function is of Heitler-London form constructed out of the atomic KS or- 
bitals [14, 15]. The physical nature of this peak structure and its connection to the Hubbard 
correlations is the main subject of our paper. We prove that the whole spatial dependence 
of the KS potential in the strongly correlated dissociation limit, including both the peak 
structure and the step (for heteroatomic molecules), is universal. It depends only on the 
asymptotic behavior of the density of the fragments, which, in turn, is mainly determined by 
the atomic ionization potentials [9] . In particular, we derive an analytic formula that allows 
us to recover the exact form of the xc potential in the dissociation limit from the knowledge 
of the ionization potentials of the independent fragments. This result adds one more item to 
the list of exact properties of the KS system and xc potential, such as Koopman's theorem, 
the exact asymptotic form of Vxc for finite systems, and the exact relation of the asymptotics 
of the density to the asymptotics of the highest occupied KS state [9] . We also demonstrate 
that the peak structure in v^^c can be viewed as a manifestation of the Hubbard on-site cor- 
relations at the level of noninteracting KS particles. The physical significance of this peak 
is that it suppresses the quantum tunneling of KS particles between two fragments, exactly 
what the Hubbard repulsion does for real electrons. This ensures that the fragments become 
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physically independent. We emphasize that this effect is not accounted for by any of the 
currently available functionals and constitutes a stringent test for the future development of 
static and time- dependent functionals aimed at describing strongly correlated systems. 

The structure of the paper is as follows. In Sec. II we discuss the physics of the strongly 
correlated dissociation limit in terms of both Hubbard on-site correlations and the KS formu- 
lation of DFT. Using a simple analytically solvable model for a 1-dimensional (ID) symmetric 
diatomic we derive the asymptotic form of the KS potential and verify our findings numeri- 
cally for more general symmetric ID systems. In Sec. Ill we uncover the universal physics 
that governs the behavior of the KS potential in the dissociation limit, derive general exact 
analytic formulas valid for all two-electron systems and verify them numerically for model 
ID heteroatomic molecules. We also discuss generalizations of the results for more general 
many-electron systems. We then conclude the paper by summarizing our main results. 

II. TOWARDS DFT IN THE DISSOCIATION LIMIT 

A. Physics of the dissociation limit: Real electrons vs. Kohn-Sham particles 

Let us first consider the qualitative physics of the dissociation of simple diatomic 
molecules with a single cr-bond formed by a pair of electrons originating from the atomic 
valence orbitals. Specific examples for this scenario include H2, Li2, and LiH, to name a 
few. When the molecule is stretched the gain in the kinetic energy due to the delocalization 
of the electrons, which is proportional to the hopping matrix element t, decreases exponen- 
tially. On the other hand, the loss in the interaction energy, due to the presence of two 
electrons on the same atom, saturates at a certain value of the Hubbard on-site repulsion 
U. Starting from some distance, roughly determined by the condition U ^ t, the on-cite 
Coulomb correlations block the inter-atomic tunneling, the electrons get localized on their 
own atoms, and the molecule dissociates into two physically independent fragments. 

Within DFT the real interacting system is modeled by an artificial non-interacting KS 
system with the same ground-state density. The non-interacting particles are subject to an 
effective potential via the KS [2] equation (atomic units are used throughout the paper) 



Since the KS particles are noninteracting there is no way to localize them on a particular 
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atom, independent of the distance d between the fragments. For a symmetric molecule, 
like H2 or Li2, the KS particles responsible for the formation of the bond always occupy a 
symmetric orbital with a probability of 1/2 to find either particle on each atom. 

Apparently, the behavior of the KS particles is very different from that of real physical 
electrons. The difference between the real word and an artificial word of KS particles be- 
comes especially striking in the regime of strong correlations, and the dissociation of simple 
molecules provides us with a bright example of this phenomenon. However, a certain phys- 
ical information, namely the ground-state density, is reproduced exactly by the KS system. 
Therefore, the real physics should be reflected in the properties of the KS system. Estab- 
lishing a map of the physics governed by the strong Hubbard on-site correlations to the 
properties of the KS potential, i.e. the map of the real word to the world of KS particles, 
is the main subject of this work. In order to find this map we mainly concentrate on a 
minimal model that captures all key physics of dissociation - the system of two electrons 
in a potential formed by two nuclei/potential wells. In Sec. HI we argue that the main 
conclusions are transferable to a more general many-electron case. 

In the case of two electrons in a singlet state only one spatial KS orbital is occupied. 
Therefore, the density is given as n{r) = 2|(y9i(r)|^ = 2ipl{r), because the orbital can always 
be chosen to be real. Hence, from inverting Eq. (1), the exact KS potential is given by 



with n{r) being, by construction, the exact ground- state density of the two-electron system. 
Hence, given the exact two-body ground-state wave function \&(ri, one can calculate the 
density n{r) = J (ir2|\l/(r, r2)p, and then recover the exact KS potential by inserting n(r) 
into Eq. (2). This formally maps the physical two-body wave function to the KS potential. 
However, extracting the physics behind this formal map is not as simple as one may think 
since in a general 3D case the wave function \l/(ri,r2) is a complicated object given fully 
numerically, e.g. from CI calculations, and, moreover, may be numerically problematic for 
realistic systems when one reaches the dissociation limit . Therefore, it is instructive to look 
first at some simplified models and then, after the essential physics is understood, return to 
reahstic situations. 

An obvious simplification, which still contains all physical ingredients of the original prob- 
lem, is to consider a system of two interacting particles in one dimension. The corresponding 
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two-electron Schrodinger equation takes the form 

+ ^^cxt(a;i) +Vc^t{x2) +Vint{\Xl -X2\) 









H 




dxl) 



'^{xuX2) = E^{xuX2), (3) 

where Vext{x) is the external potential, and fintd^; — x'\) is the potential of the inter-particle 
interaction. At the end of this section and in Sec. Ill we present the results based on the full 
numerical solution of Eq. (3). However, to gain some physical insight into the shape of Vg 
in the dissociation limit we simplify the model even further to make it analytically solvable. 



B. Analytical model of strongly correlated electrons 

First, we assume that the external potential in Eq. (3) is given by a sum of two 5-function 
wells of equal strength, v, located at the points x = ±d/2. Similarly, we take the interaction 
to be a zero-range delta-potential of strength A 

^ext(x) = -v[5{x - d/2) + 5{x + d/2)], (4) 
Vint{\x-x'\) = X6{x-x'). (5) 

Physically, in the dissociation limit the only role of the interaction is to block the inter-atomic 
tunneling. Therefore, in that limit, the behavior is expected to be universal and independent 
of a particular form and/or strength of the interaction. This leads us to the last simplifying 
assumption, namely the limit of infinitely strong 5-repulsion, A — oo. Now the problem 
becomes immediately solvable by the so called Girardeau mapping [16] (see also a more 
recent review [17]), which allows to map the ground state of strongly interacting "hard- 
core" bosons (a symmetric wave function) to the ground state of noninteracting fermions 
(antisymmetric wave function). In our two-particle case the exact ground-state (singlet, i.e. 
symmetric) wave function takes the form 

^(xi,a;2) = |0i(xi)</)2(x2) - 02(a;i)0i(a;2)|, (6) 

where (f)i{x) and (p2{x) are the two lowest states of the following one-particle Schrodinger 
equation 

-^(Pni^) + Vc^t{x) (j)n{x) = e„0„(a;). (7) 

In other words, the ground state of two infinitely interacting particles in a singlet state 
is given by the modulus of the ground state wave function of two noninteracting spinless 
fermions in the bare external potential Ve^t{x). 
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The two lowest energy solutions of Eq. (7) with the external potential of Eq. (4) are easily 
found to be (pi{x) = (p+{x) and 4>2{x) = <t>~{x) with 

M^) = C± (e-"±l^+'^/2| ^ g-a±|x-d/2|^) ^ 

where C± are the normalization constants. The parameters a±, which determine the cor- 
responding eigenvalues ei,2 = e± = — tt±/2, are the solutions of the following dispersion 
equations [22] 

a± = v{l± e~°±'^) . (9) 

Using the exact ground-state wave function (6) we obtain the exact density 

nW^/..V(...')^,^W+,iW. (10) 

and, finally, by inserting n{x) into the ID version of Eq. (2), the exact KS potential for our 
strongly correlated two-particle system, Vs{x) = Avs{x) + Vc:>,t{x), 

Equation (11) gives the exact KS potential for any distance between the wells. In the 
dissociation limit, vd ^ 1, a± ^ v and e± —v'^/2. Therefore, the last two terms in 
Eq. (11) cancel while the remaining first term simplifies to 

Avs(x) = ^- = ^ (12) 

2cosh2(2t;x) cosh2(2v^x) 

I = is the ionization potential of a separate fragment, the delta-potential of strength 

V. Hence, we have found that the exact KS potential in the dissociation limit has the form 

of a "wall" built up between the two fragments of the "molecule". The shape of this wall 

looks quite close to the peak structure observed numerically in previous works [10-13]. 

C. ID model for homoatomic dissociation 

It is physically plausible to expect that the behavior in the dissociation limit is indepen- 
dent of the particular form and strength of the interaction, and that the asymptotic form 
of Vs{x) for more general systems is similar to that given by the simple formula (12). We 
now verify this expectation for a ID system of two particles in a more general, but still 
symmetric, external potential, namely 



-'ext 



X 



1 1 

+ 



cosh^(a;-rf/2) cosh^(x + d/2) 



(13) 



The two particles are interacting with a finite range interaction potential of the form 

Vinti\x-y\) = — -. (14) 

cosh [X — y) 

The choice of the 1 / cosh^ shape of the wells and the interaction potential is arbitrary. It is 
simply a matter of convenience as the ID Schrodinger equation with a 1/ cosh^ potential is 
exactly solvable [19], which allows us to control the accuracy of our numerical calculations. 
In addition, the finite-range interaction (14) allows us to reach the dissociation limit in a 
controllable way without numerical instabilities. 

For the numerical solution of Eq. (3) with general fext and fint, we note that the ID 
two-particle problem defined by Eq. (3) can be formally interpreted as a 2D one-particle 
problem with the Hamiltonian 
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+ 



+ v!^,{x,y), (15) 



where the effective 2D one-particle potential is defined as 

""ext (^' y) = ^cxt (x) + l^ext (v) + ^'int ( 1 2^ - ^ | ) • (16) 

Consequently, the exact ground-state wave function \l/(x, y) and the exact one-dimensional 
ground-state density for the physical two-particle system, n{x) = J dy\'${x,y)\'^, can be 
obtained numerically from any computer code that is able to treat non-interacting electrons 
in two dimensions. All our calculations in this work were carried out with the OCTOPUS code 
[20]. 

The exact KS potential, Vs{x), for Vcxt and fjnt of Eqs. (13) and (14) with v = 0.9 and 
b = 0.5, and varying interwell distance d is shown in Fig. 1. At first sight the results look 
very surprising: starting from a certain distance, d = 8 a.u. for these particular parameters, 
the shape of the KS potential saturates exactly at the form given by the analytic formula 
(12) with / being the ionization potential of a single 1/ cosh^ well. Calculations for different 
strengths of the wells v, different interaction strength b, as well as for a long-range soft- 
Coulomb inter-particle interaction all show the same result [23]. At large distances the 
exact KS potential is not only similar to the analytic form of Eq. (12), obtained from 
an oversimplified model with an infinite (5-repulsion, but matches it exactly as soon the 
dissociation limit is reached! In the next section we discuss a deep, though simple physical 
reason for this seemingly surprising universality. 
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FIG. 1: Kohn-Sham potential for two equivalent wells at different distances {v = 0.9). The external 
potential has been subtracted to facilitate the comparison. Analytical results are given by Eq. (12). 

III. EXACT KOHN-SHAM POTENTIAL IN THE DISSOCIATION LIMIT 
A. Universality of the Kohn-Sham potential 

In order to understand the nature of the universal peak in the asymptotic form of the 
KS potential we turn back to our first simple model with an infinite zero-range repulsion 
and look more closely at the behavior of the exact density determined by Eq. (10). In the 
dissociation limit, vd ^ 1, the functions </'+(x) and 4>-{x) become simple symmetric and 
antisymmetric combinations of "atomic" orbitals. Taking the squares and summing them 
up, as suggested by Eq. (10), we find that all interference terms, i.e. the cross-product 
of different atomic orbitals, cancel, and the total density reduces to a sum of two atomic 
densities 



This is exactly what Hubbard on-site correlations do - they destroy the inter-atomic tun- 
neling/interference, which localizes the electrons on separate sites, and eventually makes 
the density to be the sum of the densities of two physically independent fragments. On 
the KS side of the mirror, the KS potential, whatever it is, cannot localize the KS parti- 
cles. However, by building up a self-consistent wall between the fragments it suppresses 
the tunneling/interference of the atomic KS orbitals to mimic the density distribution of 
the two independent atoms. Thus, the physics of the on-site Hubbard correlations in the 
real world is mapped to the wall in the KS potential in the artificial world of KS particles. 
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It is, therefore, not surprising that the universahty of the physics in the dissociation hmit 
is reflected in the universal form of the asymptotic KS potential. Since for DFT only the 
density distribution is essential, the general condition that determines the KS potential in 
the dissociation limit is simply 

n{r) = ni{r) + n2{r) , (18) 

in other words, the total density n{r) is equal to the plain sum of the densities, ni{r) and 
n2(r), of the two independent fragments. The asymptotic form of the KS potential should 
be such that it supports the density distribution given by Eq. (18). As the densities ni(r) 
and n2(r) decay exponentially from different sides the only way to mimic this at the level of 
a single KS orbital is to insert a potential wall in the middle region. 

Having understood the key physics we are ready to go to more complex systems. 



B. Kohn-Sham potential of heteroatomic ID molecules 

It is now straightforward to find the form of the KS potential in the dissociation limit for 
a general ID molecule formed by two different wells. Assuming that the densities, ni{x) and 
n2{x), corresponding to one electron sitting in a separate well are known, we require that 
the total density n{x) is given by their sum, Eq. (18), and substitute this sum into Eq. (2). 
The result can be reduced to a form that looks structurally similar to Eq. (11) 

2(m + n2)2 + n,+n2 ' ^ ' 

where /1.2 are the ionization potentials of the fragments and / = min{/i, I2} is the ionization 
potential of the total system. Equation (19) is valid in the dissociation limit, and from its 
structure it is clear that Avs{x) has a nontrivial x-dependence (i.e. differs from a constant) 
only far away from the "atoms" , where the densities fall off exponentially. Therefore, for the 
practical evaluation of Vg in the dissociation limit, it is sufficient to know only the asymptotic 
behavior of the density of the separate fragments. In the ID case, the asymptotics of the 
densities ni{x) and n2{x) have the following general form 

ni,2(x) = Ai,2e-'"^'^l"±'^/'l, (20) 

where the exponents ai^2 are related to the ionization potentials of the atoms Ii^2 = (^12/'^ 
and y4i 2 are prefactors to the exponential decay. 
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Inserting Eq. (20) for a symmetric molecule (equivalent wells with Ai = A2, and ai = 0:2) 
into Eq. (19) we immediately recover our first model result of Eq. (12) thus confirming its 
universality. In a general asymmetric case (different wells or a "heteroatomic" molecule) 
a new qualitative feature, a "shelf", appears. This shelf in Vg is such that it aligns the 
ionization potentials of the atoms. Formally, it results from the last two terms in Eq. (19) 
which do not cancel if the ionization potentials are different. Substituting the general form 
of Eq. (20) into Eq. (19) we find the explicit results for Vg in different regions of space. 
It is convenient to represent the KS potential as a sum of two different contributions, i.e. 
Avs{x) = vi^\x) + vi!^\x). For the region between the wells, i.e. for —d/2 < x < the 
two contributions correspond to the "wall" and the "shelf" discussed before. They are given 
by 

cosh^[(ai + 02) (x + xo)] 



-^\-) = __.2r . - (21) 



^f\^) = , ^ (22) 

l + exp[2(ai + a2)(x + a;o)] ^ ^ 



with 

1 

Xq 



ai + ^2 



"1 - «2 , , , A2 



(23) 



Here, and in the following, we have assumed that ai < 02, i.e. that the left fragment 

(2) 

has a larger ionization potential. Obviously, for a symmetric configuration, Vg vanishes 
identically, i.e. there is only a peak in this case. Also, in this case xq = 0, i.e. the peak is 
exactly in the middle between the two identical fragments as expected from symmetry. 
For X < —d/2 the two contributions read 

vi'\x) = , ~ , (24) 

cosh [(«! — 0:2) (2; + Xq)] 

' ^ ' l + exp[-2(ai-a2)(x + x^)] ^ ' 



with 

' '"^ + "^.rf-los^l. (26) 



Xn — 



ttl — 0^2 



2 ° Ai 



Contrary to before, v^s^ does not describe a peak but it can actually be shown that the 
potential is strictly monotonically increasing describing the building up of the shelf, or its 
return to zero depending on the direction one approaches x'q from. Also, for the symmetric 
case, both contributions vanish as there is no shelf in that case. For the region x > d/2 the 
potential decays exponentially without specific features. 
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We emphasize that neither the specific form of the fragments nor the type of interaction 
between the electrons enters the derivation of our analytical result directly. The specifics 
of the fragments appear in the result only via the parameters ai^2 and Ai 2. The former 
describes how fast the density decays, i.e. it is directly related to the ionization potential of 
each fragment. The latter is connected to the normalization of the wave function and only 
enters the potential as a logarithmic correction to the position of the peak and shelf without 
changing the shape of the potential. For a symmetric system the potential is completely 
determined by the ionization potential of the two fragments. In all cases, symmetric and 
asymmetric, the functional form of the KS potential is universal, only the position and 
the width and height of the peak depend on the system under consideration. Both in the 
symmetric and in the asymmetric case the presence of the universal wall reflects Hubbard 
correlations. The potential wall suppresses the tunneling and drives the KS density to the 
density corresponding to physically independent subsystems. 

To ensure that our universal analytical formulas are indeed correct we performed numer- 
ical calculations for an asymmetric two-electron system with the external potential given by 
the sum of two different potential wells 



Vl V2 

cosh^(x - d/2) cosh^(x + d/2) 
with Vl = 0.9 and V2 = 0.7. As before, for the interaction we keep the finite range potential of 
Eq. (14). Since the one-particle problem with 1/ cosh^ is exactly solvable [19] the parameters 
ai.2 and A1.2, entering our asymptotic formulas Eqs. (21)-(26), are available in the analytic 
form. In particular, for the pre-exponential factors in the asymptotics of the "atomic" 
densities we get 



.2....rK2 + i/2) 



V7rr(ai,2) 

where F denotes the usual Gamma-function. 

In Fig. 2 we show the comparison of the analytic KS potential given by Eqs. (21) and 
(22) and the KS potential obtained from the full numerical solution of the problem defined 
by Eqs. (3), (14), and (27). As expected, in the asymmetric case, Vi 7^ V2, the KS potential 
acquires a shelf structure in addition to the peak. The shelf is a direct result of the necessary 
alignment of the KS energy levels (the ionization potentials) in the two fragments [9, 21]. 
It is already not so surprising to see that the KS potential again approaches the analytic 
asymptotic form with increasing distance. As in the analytic calculation, the exact position 
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FIG. 2: Kohn-Sham potential for two different wells at different distances {vi = 0.9, V2 = 0.7). The 
external potential has been subtracted to facilitate the comparison. Analytical results are given 
by Eqs. (21) and (22). The potential returns to zero at large negative x. 

of the peak and the shelf depends slightly on the distance between the two wells always being 
closer to the deeper well. While in the symmetric case, see Fig. 1, the dissociation limit is 
reached at a distance of around 8 a.u. in the asymmetric case around 11 a.u. are necessary. 
In both cases the numerical results agree perfectly with the analytical expression. The larger 
distance, necessary in the asymmetric CcLSG, IS db result of the shallower right potential well 
in that case. 

Unfortunately, the analytical result of the shelf returning to zero can not be verified 
numerically for the systems at hand. The position Xq, Eq. (26), is so far away from the 
actual potential wells that the density is numerically zero. There is, however, no doubt that 
the shelf returns to zero exactly as predicted by the analytic formula. 

C. Generalizations to three-dimensional and many-electron systems 

The general argumentation used in the previous subsection to derive the exact KS poten- 
tial in the dissociation limit is not restricted to ID systems. The general physical condition 
for dissociation is that the density is given by the sum of the densities of the independent 
fragments, Eq. (18), because the inter-fragment tunneling is destroyed by Coulomb correla- 
tions. The inversion formula of Eq. (2) is also valid for any two-particle system independently 
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FIG. 3: Kohn-Sham potential of the hydrogen molecule in the dissociation limit. 



of dimension. Therefore, an elementary 3D generalization of Eq. (19) takes the form 

where all notations are the same as in Eq. (19). Using this formula we can recover the 
exact limiting functional form of the KS potential for any two-particle object dissociating 
into two one-particle fragments. The only required input is the long-range asymptotics 
of the independent fragments, which is mainly determined by their ionization potentials. 
It is important to emphasize that the pre-exponential factors give only weak logarithmic 
corrections to the position of the wall and the shelf. 

As an illustration, we present the exact KS potential that controls the dissociation limit 
of the H2 molecule. The final results obtained by inserting the ionization potential of the 
hydrogen atom, and the electronic densities of two independent hydrogen atoms, located at 
the points Ri and R2, into Eq. (29) takes the form 

A.f»(r) = (30) 

4 cosh (ri — r2) 

where ri^2 = r — R1.2 are the vectors between the two protons and the considered point in 
space. The KS potential for the hydrogen molecule, Eq. (30), is shown on Fig. 3. It is easy 
to see from Eq. (30) that along the molecular axis Av^^{p = 0, z) is exactly of the ID form 
Eq. (12), while in the perpendicular direction it has a Lorentzian shape, 



Avf'ip, z = 0) = 1 , (31) 



with the width increasing at increasing distance between the two hydrogen atoms. 

Similarly, we can obtain an explicit form of the exact KS potential for any two-electron 
system in the strongly correlated dissociation limit. Moreover, one can argue that the general 
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formula (29) remains valid also for many electron systems in those cases where the separate 
fragments have a single electron in the highest occupied KS orbital. Indeed, in this case the 
asymptotic behavior of the density away from the atoms is completely determined by the 
two KS particles in the highest occupied KS molecular orbital (KS HOMO), while the rest 
of the electrons effectively contribute to the rigid atomic cores. Therefore, the asymptotic 
form of the KS potential can be obtained by inverting only one KS equation, namely for the 
KS HOMO, and, hence, the two-particle formula (2) remains asymptotically valid. 

IV. CONCLUSIONS 

In conclusion, we have presented a recipe to calculate the exact KS potential of systems 
in their dissociation limit. The main ingredient is the ionization potential of the dissociated 
fragments, a quantity that is readily available from spectroscopic data. We have presented 
the explicit results for a one-dimensional model system and the hydrogen molecule. It is 
shown that the functional form of the potential is independent of the specific system and 
the details of the interaction as long as the latter is repulsive and sufficiently strong. For the 
ID model system the numerical results approach the analytical one as the distance between 
the two fragments is increased. Hence, they confirm our analytical result perfectly for both 
a symmetric and an asymmetric system. Our results not only pose a strong constraint for 
the development of exchange-correlation functionals but also introduce an alternative way 
to look at the electron localization in strongly correlated systems. How to incorporate those 
effects in a density-functional treatment remains a challenge. It is especially intriguing to 
explore implications of our universal results for the quantum transport in the regime of 
Coulomb blockade. It is natural to expect that the potential wall in the KS potential should 
modify the tunneling probability when the transport is described in terms of KS DFT. 
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